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We consider a phase-space model for particle production in nuclear collisions. Once the mul- 
tiplicities of the individual particle species are known, single-inclusive momentum spectra can be 
computed after making simplifying assumptions for the matrix element for multiparticle production. 
Comparison of the calculated spectra with data for pions and kaons from central Pb-|-Pb collisions 
at -ELab = 158 AGeV reveals a residual longitudinal phase-space dominance in the final state of 
the reaction. We account for this by modifying the isotropic, relativistic invariant phase space in 
a way which retains boost invariance in beam direction but suppresses large transverse momenta. 

. . . Adjusting a single parameter, we obtain a reasonably good description of transverse momentum and 

lO rapidity spectra for both pions and kaons. 
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O . I. INTRODUCTION 

o 

\^ • In relativistic heavy-ion collisions one can study strongly interacting matter under extreme conditions [1]. Of 
^sj ' particular interest is the question whether a novel phase of nuclear matter, the so-called quark-gluon plasma (QGP), 
can be created in such collisions. In this phase, quarks and gluons are deconfined and chiral symmetry is restored. 
CO ' Lattice QCD calculations [2] predict that, in thermodynamical equilibrium, this phase occurs at temperatures above 
^ 170 MeV (for vanishing baryon chemical potential). Such temperatures are expected to be reached and exceeded 

for collisions at CERN-SPS and BNL-RHIC energies, respectively. 
, In nucleus-nucleus collisions at ultrarelativistic energies, secondary particles are copiously produced over a small 
' region in space-time. It is then natural to assume that interactions between these particles drive the system towards 
' (local) thermodynamical equilibrium. Indeed, measured single-inclusive transverse momentum spectra of individual 
• particle species show an exponential decrease similar to that expected in thermal equilibrium. Following this line 
' of arguments, one has studied the evolution of the system using models based on ideal fluid dynamics where the 
system is assumed to be in local thermodynamical equilibrium. At RHIC energies, such models are able to reproduce 
the data for single- inclusive transverse momentum spectra of various particle species [3] . A particular success of the 
fluiddynamical approach was the correct prediction of the magnitude of elliptic flow. The fact that data on elliptic 
flow seem to saturate the ideal fluid limit has recently led to the speculation that matter created in collisions at RHIC 
(— I ] energies is a so-called strongly coupled QGP with very small viscosity over entropy ratio [4]. 

However, the fluiddynamical approach is not without problems. A long-standing problem is the so-called "HBT puz- 
zle", i.e., the fact that HBT correlation radii and single-inclusive particle spectra cannot be simultaneously described 
within the fluiddynamical approach [5]. 
' Another problem of the fluiddynamical approach concerns the reproduction of the measured particle ratios. Under 
the assumption of local thermodynamical equilibrium, as well as conservation of electric charge and strangeness, one 
can determine from such particle ratios the temperature and the baryon chemical potential of the system on the 
hypersurface where chemical freeze-out occurs, i.e., where inelastic collisions cease and the multiplicities of produced 
particles no longer change [6-8]. The values for the chemical freeze-out temperature and baryon chemical potential 
thus obtained indicate that the system undergoes chemical freeze-out well before thermal freeze-out, i.e., when also 
elastic collisions cease and the single-inclusive particle spectra no longer change their shape. If one accounts for this 
in the fluiddynamical evolution by allowing for departures from chemical equilibrium [9] and thus keeps the particle 
ratios at their experimental values, one is no longer able to describe data for collective flow [10]. (An exception is the 
study of Ref. [11] which reproduces the transverse mass spectra and multiplicities for pions, kaons, and (anti-)protons 
for collisions at RHIC energies assuming simultaneous thermal and chemical freeze-out.) 

In this paper, we pursue a different strategy of describing single-inclusive momentum spectra of particles produced 
in heavy-ion collisions. Our arguments are solely based on the phase space of the particles in the final state. Thus, 
we follow a well-known and well-documented approach in the physics of hadron-hadron collisions [12-14] and apply 
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it to nucleus- nucleus collisions [15 -19]. In order to determine the phase space, we need to know the experimentally 
measured multiplicities of secondary particles. We also have to make simplifying assumption for the matrix element 
of multiparticle production. We can then determine single-inclusive particle spectra and compare them with the 
experimentally measured spectra. 

This phase-space model to describe single-inclusive particle spectra represents an approach with the minimum 
number of assumptions concerning the dynamics of the reaction. In its most simple form it relies solely on the 
constraints imposed by the kinematics of the reaction, such as total energy and particle multiplicities. Deviations 
from the phase-space model are a clear indication for nontrivial dynamics beyond pure kinematics. 

Our approach is complementary to the fluid- / thermodynamical approach. There, the exponential decay of transverse 
mass spectra as seen in the data is the motivation to argue that matter is in local thermodynamical equilibrium. On 
the other hand, the data deviates from this exponential decay, especially in the region of small transverse momenta. 
These deviations are commonly attributed to resonance decays and to collective motion. Particle ratios can be 
predicted, once temperature and baryon chemical potential at chemical freeze-out are known. On the other hand, 
in our approach particle multiplicities serve as input to determine the size of multiparticle phase space. Once they 
are known, single-inclusive particle spectra can be computed. As we shall see, the calculated transverse mass spectra 
also decay exponentially. Similar to the fluid-/thermodynamical approach we then interpret the deviation of the data 
from a pure exponential decay as a signature for nontrivial dynamics, such as resonance decays and collective flow of 
matter. 

For this study, as an example we use data [20-25] of the NA49 collaboration for the 5% most central Pb-I-Pb 
collisions at a beam energy of E'Lab = 158 AGeV, corresponding to y^iiVTV = 17.3 AGeV in the center-of-momentum 
frame. For the 5% most central collisions, the average number of wounded nucleons is (iV^) — 362. Therefore, the 
total energy deposited in the reaction is estimated to be y/s = ^s/vw-^^^ = 3131.3 GeV. In principle, we could also 
have used data from other experiments, but the advantage of the NA49 experiment is that it provides full rapidity 
coverage. It thus allows us to investigate rapidity spectra in addition to the transverse mass spectra. 

This paper is organized as follows. In Sec. II we define and compute the phase-space integral. In practice, the com- 
putation is not possible unless one makes simplifying assumptions for the matrix element for multiparticle production. 
We investigate a product ansatz for the matrix element, which allows to factorize the phase-space integral and thus 
facilitates the computation. We then study two special cases in more detail, one which is completely isotropic in 
phase space and preserves general Lorentz invariance, and another which incorporates the longitudinal phase-space 
dominance of the reaction and only preserves Lorentz invariance under boosts in beam direction. In Sec. Ill we 
compute single-inclusive momentum spectra, i.e., transverse mass and rapidity spectra, from the phase-space integral 
and point out differences to momentum spectra in thermal equilibrium. In Sec. IV we confront the calculated spectra 
with experimental data. To this end, it was also necessary to reconstruct the multiplicities of unmeasured particle 
species using symmetry considerations, an empirical rule, and a simple quark model. Section V concludes this work 
with a summary and an outlook. Our units are h = c = ks = ^ and the metric tensor is chosen as g^^ = (+,—,—,—). 



II. THE PHASE-SPACE INTEGRAL 



Let us consider the cross section for production of secondaries in collisions of elementary particles (or nuclei). The 
particles (or nuclei) in the incoming channel have masses nia, rrib and four- moment a Pa, pi,. The centcr-of-momentum 
energy (squared) of the reaction is s = P^, where = {E,!*) = Pa + Pb is the total four-momentum. We assume 
that there are n particles in the final state with masses mi, . . . m„ and four- momenta Pi,- ■ - Pn, where = {Ei, p,), 
with Ei = y^pf + m? . The cross section a for this reaction is given by 

CT=^In{Pn, (1) 

where 

F = 2 ^/X{s,ma,mb) (27r)^"-^ (2) 

is the flux factor, with A(.t, y, z) = + y'^ + — 2xy — 2yz — 2zx being the kinematical function [13] . /„ is an integral 
over the 3n-dimensional invariant momentum space of the secondaries, 

inipn = [fl^ ^'(p" Ep") nPi) ■ (3) 

The integrand consists of a four- dimensional, four-momentum conserving ^-function and of the square T{pi) of the 
transition matrix element. 
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The secondaries produced in a reaction come in different species r. The multiplicity of particles of species r is 
denoted by n^, such that n = ^^n^. In order to compute /„ one needs to specify the multiplicities n,.. These differ 
from event to event. Usually, event-averaged quantities arc given in the literature. We shall exclusively use these 
event-averaged multiplicities in the following. We shall see that our results do not change qualitatively if we vary the 
multiplicities Ur within a factor of where we assume to be a measure for the statistical error. One could 

also study multiparticle production on an event-by-event basis [26] in order to examine particle number fluctuations, 
but this is beyond the scope of the present work. 

In order to calculate /„, according to Eq. (3) 3n integrations have to be performed, where in our case n is larger than 
2000. The (5'*-function prevents that these integrals decouple. Therefore, the computation of /„ requires in general 
numerical methods. One can distinguish three approaches, recursion relations [12, 13, 27], Monte Carlo methods 
[28, 29], and the statistical approach [13]. In this work we employ the latter, because they are most appropriate for 
a large number of secondaries. The essential idea is to use standard methods of statistical mechanics [30-35]. In 
statistical mechanics, calculations performed in the microcanonical ensemble are difficult, because the energy is fixed. 
Therefore, one replaces the microcanonical ensemble by the canonical ensemble, where the energy can vary, but the 
temperature is fixed. In the thermodynamical limit, i.e., for large particle numbers, both approaches yield identical 
results. The same idea is applied here to compute /„, but there exist several methods which differ in detail. Some use 
the central-limit theorem of probability theory, as suggested and explored by Khinchin [36] . In this work, we employ 
the saddle-point method [34], which yields, in the approximation used here, the same result as the methods based on 
the central-limit theorem. 



A. The calculation of the phase-space integral 

In general, for multiparticle reactions the matrix element T(pj) is too complicated to calculate analytically. There- 
fore, in order to make progress one resorts to simple approximations for T{pi). Such an approximation can be 
motivated by the symmetries of the problem, and it should not spoil the Lorentz invariance of /„. Our particular 
ansatz used in the following is 

n 

T{Pi) = l[fi{p±i). (4) 

The most simple case, without any assumptions about the symmetries of the problem and preserving full Lorentz 
invariance, is 

fi{p±i) = 1 • (5) 

Another case studied here is motivated by the longitudinal phase-space dominance in hadron-hadron and high-energy 
nuclear collisions. Choosing 

fi{p±i) = exp [-a{m±i - m,)] , (6) 

where m±i = \Jv\i + is the transverse mass, the transverse phase space is exponentially suppressed. This ansatz 
preserves boost invariance with respect to the beam axis, but is no longer invariant under Lorentz transformations 
perpendicular to this direction. 

Inserting the ansatz (4) into Eq. (3), we obtain 

= / n ^ /^(p^o ^\p' - ") • (7) 

The four-momentum conserving delta-function introduces a mutual dependence of the various momentum integrals 
and thus makes the integral /„ difficult to compute directly. This problem is similar to that which one encounters 
in the computation of the microcanonical partition function in statistical mechanics. In that case, the computation 
is facilitated by considering the canonical partition function as the Laplace transform of the microcanonical one. 
Technically, the Laplace transformation decouples the momentum integrals. Once the canonical partition function 
is known, an inverse Laplace transformation brings us back to microcanonical partition function which we wanted 
to compute originally. In our case, the analogue of the canonical partition function in statistical mechanics is the 
(four-dimensional) Laplace transform of the integral /„, 



= / d^P exp(-a,P^) 1^{P^) 



(8) 
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The Laplace transform is a function of the four- vector a^. Inserting /„ into Eq. (8) and performing the integration 
over P'^, we obtain 



= /nS^/i(P^Oexp(-a.pr)- (9) 

i=l ' 



We define the functions 



such that 



2^ fi{p±i) exp (-a^pH , (10) 



^n{a^)=llM'^f.) ■ (11) 
The original integral 7„ can be obtained from the inverse Laplace transformation, 



1 pc^+ioo 

^"(^'') = TT^i / exp (a.P^^) $„(a^) . 



(12) 



The path of integration is parallel to the imaginary axis, where is a constant 4-vector. Each component is chosen 

such that the path of integration lies in a range of values for . where the integrand does not have any singularities. 
We proceed by taking into the argument of the exponential function, 

^"^^"^^(2^/ . ^'"exp[F(a^)] , (13) 

where we defined 

F{a^) = a,P'' + \n<^,,{a^) . (14) 

In general, the integral (13) cannot be solved analytically. Statistical methods are used to obtain an approximative 
solution. Here, we employ the saddle-point approximation [13, 34]. The saddle point of the integrand, is given by 
the condition 



= , (15) 



dai, 

or, with Eq. (14), 

ain$„(/3^) 1 /■ d^Pt ft ^ 7. / A\ 

Note that, in the statistical mechanical problem of computing the microcanonical partition fmiction, the parameter 
Pn would correspond to the inverse temperature. Going into the rest frame of the saddle point, = (/3, 0) and using 
the fact that f{p±i) and p^ = Ei are even functions of momentum, we observe that this rest frame is identical to 
the centcr-of- momentum frame, = (-^5,0). Thus, and P^ must be proportional to each other, P^^ ^ /?^, in 
any frame which differs from the rest frame by a boost in longitudinal direction. We now determine the constant of 
proportionality. To this end, we compute the four-gradient of In at the saddle point assuming to be a function 
of the absolute value /? = a/A^ of Z?**, and its direction, = f)^'/|3. Then, 

ain$„(/?,/3^) a In $„(/?, /3^) ain^„(/3,/3^) 1 ^ , 

= — w. — = — m — ^ — — 



where 



Bx'' = l3^=gx''-M'' (18) 
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is the projector onto the direction which is four-transverse to However, since f3^ and P^^ arc proportional to each 
other, the part proportional to Bx" must be zero, i.e., there is no dependence of ln$„ on /3^, and we have 



gin 

d(3 



Going to the center-of-momentum frame, this equation serves to (implicitly) determine the absolute value of (3, 

51n$„(/3) _ 1^ 1 f,s„ ,^.f,E. 



(19) 



(20) 



This equation relates (3 to the center-of-momentum energy ^/s. In the analogous problem of statistical mechanics, such 
a saddle-point condition relates the inverse temperature (the independent thermodynamical variable in the canonical 
ensemble) to the total energy (the independent thermodynamical variable in the microcanonical ensemble). 

We now proceed with the evaluation of the integral (13). We choose the integration path to lead through the saddle 
point and change the integration variable to = —i{a^^ — /?^). For large particle numbers the integrand is sharply 
peaked around the saddle point, = Q, and we may employ a Taylor expansion of the function F{a^). For our 
numerical calculations presented below, we checked that terms higher than second order in this expansion are at most 
of the order of ~ 10"** and thus insignificant. We therefore restrict the following consideration to terms of second 
order in t^, 



F{a^)=F{P^)--f't^F,,x 



(21) 



where we used the saddle-point condition (15) to eliminate the linear term in t^, and where the second derivative of 
F{c(fi) with respect to at the saddle point is denoted as 



Using Eq. (14) and the fact that ln$„ only depends on the absolute value of we can further evaluate F^x, 

92 In $„(/?) 



where we employed Eqs. (18) and (20) and defined 

,2 92 In $„(/?) «2 



P 



(22) 



(23) 



(24) 



For the explicit expression on the right-hand side we have chosen the rest frame of (i'^. Inserting the Taylor expansion 
into Eq. (13) we obtain 



ycitexpj-^ 
We evaluate the t-integral in the rest frame of Z?**, where 



(27r)^ 



i.e., the integral decouples into four Gaussian integrals. The final result is 



IniP^ 



a ( a 



47r2 



■$„(/3)^ \^ 



3/2 



(25) 
(26) 

(27) 
(28) 



For large multiplicities n, the logarithm of this expression can be estimated as 

In /„(P'')~/J^/i + In $„(/?) . 

Here, we have only kept terms of order 0(n), such as y^. Terms proportional to ln/3 are of order 0(1), because jS 
itself is ~ 0(1), see e.g. Eq. (34) below for the massless case. Also, terms of order In-^s and Iuk are of order 0(1), 
because ^/s and k are of order 0{n), see e.g. Eq. (35) below, and thus Iny^ ^ In/t 0(1). Thus, with Eq. (20) we 
obtain 



ain/„(P^) _ ^ ain/„(P^) dp 

= 2\/S — p + 



ds 



In $„(/?) 
8(3 



(3. 



(29) 
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B. Special cases 

We now discuss some special cases of the result (27). In the case (5), is invariant under any Lorentz transfor- 
mation and thus can only depend on the Lorcntz-scalar quantity /? = ^J(3'^l3v It may be computed in any frame, but 
the rest frame of (3^ is particularly convenient. In this frame [37, 38], 

/■°° r? 2irm ■ 

ct>m = 2-K dpi^^e-l"^^ = ^Kr{pmi) , (30) 

and $(/?) is readily obtained from Eq. (11). The saddle-point condition (20) reads 

_ 2n KoiBnij) 



and K is determined as 

2 ' 



If all particles are massless, we obtain 



(32) 



= S , (33) 



where the saddle point is given by 



Finally, 



such that 



The exact result is [34] 



/3=^- (34) 
K = An- I3s/s = 2n , (35) 



2n-4 

/s ■ 



= it) ■ w 



^■■'^'•>° 2.-Hn-l)!(n-2)! -^'""- 

For large n, one may approximate the factorials with the help of Stirling's formula. Using (1 — x/n)" for 
n ^ oc, one then observes that the result (36) obtained via the saddle-point method is consistent with the exact 
result (37). 

In the case (6), is invariant under Lorentz transformations along the z-axis. Again, we compute it in the rest 
frame of 0^. Defining the function 

/•OO 

Ji,m{x,y) = 27reM dz z^-^'/^ Km{z) , (38) 

J X 

we obtain 

<t)i{P) = JiflWrui, arm) ■ (39) 

The saddle-point condition (20) now reads 
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FIG. 1: The logarithm of the phase-space integral for n = 2611 massless particles (solid curve) and for massive particles as 
listed in Table II (dashed curve). Both curves are computed for o = 0. 



The coefficient k is determined as 



In order to get a feeling for the magnitude of the phase-space integral consider the hadron multiplicities given in 

Table II. There are n = 2611 secondaries with a total rest mass of ^ = 824.65 GcV. The average particle mass 
is {mi) ~ 0.316 GeV. The total available energy is = 3131.3 GcV, see introduction. For a = 0, this results in 
In ^ 0.5 • 10^^^^, or ln/„ ~ 6600. The saddle point of the t-integration is located at /? = 2.027 GcV"^ In Fig. 1 we 
show the logarithm of the phase-space integral as a function of ^/s for 2611 massless particles and for the set of 2611 
particles of Table II. For a given y/s, the result for massive particles is always below that of massless particles, since 
some of the energy is used to generate the rest mass of the particles and is no longer available for the phase-space 
volume occupied by the system. Only when ^/s oo, the massive case approaches the result for massless particles. 

The influence of the parameter a on the value of the phase-space integral is demonstrated in Fig. 2. The phase- 
space integral is again computed with the set of particles listed in Table II. The value of the integral decreases with 
increasing a, which is naturally explained by the exponential suppression of transverse phase space, cf. Eq. (6). 



In this section we demonstrate how one can compute single-inclusive particle spectra from the phase-space integral 
discussed in the previous section. For comparison, we also present the calculation of single-inclusive spectra in thermal 
equilibrium. 




(41) 



III. SINGLE-INCLUSIVE PARTICLE SPECTRA 
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A. Particle spectra from phase space 
According to Eq. (1), the differential cross section for inclusive production of, say, the nth secondary particle is 



En 



d^Pn F d^p. 

fn{P±n) frjd^p 



2F 

fn{P±n) 



i=l * j 
In-l{P^-pO. (42) 



2F 

In-i is the phase-space integral for the first n — 1 particles at a reduced center-of- momentum frame energy, 

s' = {P- pnf = s + ml-2E„^/s = s + ml-2 m_L„ cosh2/„ . (43) 
The single-inclusive particle spectrum is obtained from the differential cross section (42) via 

d^Pn ~ ^"""^ O d3p„ - ^""^ 2 /„(PM) • ^^^^ 

As discussed in the previous section, the calculation of the phase-space integrals on the right-hand side will be done 
in the rest frame of the collision. However, because there are n — 1 particles in the numerator and n particles in the 
denominator, the rest frame of /„_i will be different from that of /„. This is irrelevant if the phase-space integral is 
Lorentz invariant, such that In-i{P^ —Pn) = -^n-i(s') and In{P^) = /n(s). If the phase-space integral is only Lorentz 
invariant with respect to boosts along the 2;-axis, this difference is also irrelevant if the nth particle moves along this 
axis. If not, we assume the recoil in transverse direction to be sufRciently small such that the rest frame of In-i is 
approximately the same as that of /„. This is certainly a good approximation for large n. 

Transverse mass spectra are obtained from Eq. (44) after integrating over the transverse polar angle. 



(45) 
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An additional integration over m_L„ yields the rapidity spectrum, 

'^"^ / \ f°° ^ f I \ ^n-l{s') , . 

~ TT / dmj_n m_L„ /ri(p_L„) , . . . (46) 

"'Vn Jmn ^n\S) 

A further integration over ?/„ then gives the multipUcity of the n-th particle, which is an input parameter. This 
integration serves as a check for the consistency of our results. 

In the following, we discuss several special cases of the general result (45) for the transverse mass spectra. 



B. Special cases 

Let us consider the case (6) for the function fi{p±i). The other case (5) can be simply obtained by setting 
a = in the final result. In Eq. (45), exponentiate the numerator /„_i(s') and expand ln7„_i(s') to first order in 
s' — (s + m^) = — 2m_L„cosht/„-ys [27]. This is a good approximation, since the energy of a single particle is small 
compared to the total energy in the system. En = m_L„coshj/„ <^ ^/s. We may also drop rn^ compared to s. We 
obtain 

1 (Pnn ^, _ , , „™ /„-i(s) 



/ \ am ^n-l\S) ( ?Tlj_„ \ 



Here we defined the slope parameter 

^ = 2 \/i cosh?/„ ln/„_i(s) + o ~ /3 cosh?/„ + a , (48) 

where we have used Eq. (29) and the fact that n — 1 ~ n for large multiplicities. The transverse mass spectra exhibit 
an exponential decay in the transverse mass m±n with a slope parameter given by Eq. (48). Note that T* is not a 
free fit parameter that one could adjust to describe the data, like in global thermal fits of single-inclusive particle 
spectra. In contrast, it is unambiguously determined by the value of the saddle point /3, the rapidity y„ of the particle 
considered, and the parameter a. The slope parameter depends on the mass of the considered particle only through 
the position /? of the saddle point of the phase-space integral. However, this dependence is very weak, because the 
position of the saddle point is largely decided by the values of the total multiplicity n and the total available energy 
y/s ^ rUn ■ Quantitatively, the slope parameter varies by less than 0.1 % in the range of particle masses — 2 GeV. 
As a consequence, in this simple phase-space model the transverse mass spectra of all particle species have the same 
slope. 

In Fig. 3 we demonstrate the exponential decay of the transverse mass spectra by showing tt^ transverse mass 
spectra for different rapidities, computed from Eq. (47) with a = and using the particle multiplicities of Table II 
to compute the phase-space integrals. One observes that the slope parameter increases with increasing rapidity, as is 

obvious from Eq. (48). 

The case where a = and all particles are massless deserves a special treatment. If all particles have the same 
mass, then the mean multiplicity of the nth particle is simply the total multiplicity, (n„) = n, and we obtain for the 
transverse mass spectrum (omitting the index n) [16] 

1 <fn _^^ In-i{s') _ 2n(n-l)(n-2) / 2m^coshy y-' ^^^^ 



m_L dm±dy Inis) s \ ^/s 

where we have used Eqs. (37) and (43). Using the saddle-point condition /? = 2n/y/s, cf. Eq. (34), and (1— a;/n)" ~ 
(valid for large n), we obtain 

m_L dm±dy 2 

where we have dropped subleading terms of order 0(l/n). One observes that the transverse mass spectrum is 
exponentially decreasing as a function of the particle energy E = m± coshy. From Eq. (34) we see that the inverse 
slope parameter 1/(3 corresponds to half of the available energy per particle. We note that the result (50) could have 
also be obtained from Eqs. (47) and (48). To see this, employ a = and Eq. (37) to show that l/T* = /3. Note 
that the apparent "temperature" T* is half the energy per particle, T* = | ^/s/n [16, 27]. This is in contrast to the 
temperature T for an ultrarelativistic gas in thermal equilibrium, which appears instead of T* in the transverse mass 
spectrum, see Eq. (55) below, and which is one third of the thermal energy per particle, T = ^E/n. Thus, if we 
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FIG. 3: Transverse mass spectra for tt -mesons in Pb+Pb-reactions at 158 AGeV. 



compare transverse mass spectra for the phase-space model and the thermal model at the same total energy, ^/s = E, 
the apparent temperature is a factor 3/2 times the thermal temperature, T* = | T. The inverse slope parameter T* 
of the transverse mass spectra should thus not be interpreted as the temperature in thermodynamical equilibrium. 

Integrating either the exact result (49) or the approximate result (50) for the transverse mass spectrum of massless 
particles over m±, we obtain the rapidity spectrum 



dn 



dy 2 cosh^ y 



(51) 



The half width Ay of the rapidity spectrum for massless particles is independent of the particle number and the total 
energy ^/s of the reaction, Ay = ln(l + \/2) ~ 0.8814. A further integration gives 



h'fy=^j 



dy 
2 cosh^ y 



n , 



(52) 



as expected. 



C. The thermal model 



We shall compare our transverse mass spectrum computed from the phase-space integral with those obtained in a 
thermal model. The latter are obtained from the Cooper-Frye formula [39] 

^" — = 1 T~ = / da„p„gn{x,p„) , (53) 



d^Pn rn±n dm±ndyn Jt. 

which counts the number of on-shell particles with three-momentum p„ streaming through a space-time hypersurface 
E with normal vector dcr^. Here, gn{x,Pn) is the single-particle on-shell momentum distribution function for the nth 

particle at space-time point x and E is the hypersurface where kinetic freeze-out of particles occurs. In general, this 
hypersurface could have a complicated shape in space-time. Here, we make the simplifying assumption that freeze-out 
occurs at fixed time in the center-of-momentum frame, such that da^ = {d^x,0). Provided that gn{x,Pn) does not 
change over the freeze-out hypersurface, the integration over da,^ then simply yields a factor V, the total volume of 
the system. We also take the particles to be distributed according to the appropriate quantum statistical distribution 
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function in thermal equilibrium, _g„(a;,p„) = (i„/(27r)'^ (^p(En-ti„)/T _|_ where dn is the degeneracy and /i„ the 

chemical potential. The upper sign holds for fermions and the lower for bosons. The invariant momentum spectrum 
(53) reduces to 

1 dUn ^ dnV m_LnCOShy„ 
m±n dm±ndyn (27r)3 e('"-L»'^°s'^y"-''")/^ ± 1 ' ^ ' 

For particle energies En = m_L„coshy„ T, (U„, this becomes the classical Boltzmann distribution, 

1 dun dnVe^'-/'^ 



■ m±n COSht/„ e-™^- coshy„/T _ (55^ 



Note that this thermal distribution is not an exponential in the particle energy E„ = m_L„coshy„ because of the 
prefactor m±n cosh y„ which introduces logarithmic corrections to the pure exponential decay. 

In this work, we do not attempt to fit the temperature, the volume, and the chemical potential to reproduce the 
particle multiplicities. There is a vast amount of literature regarding such fits, see e.g. [6 8]. Our goal is to compare 
the invariant phase space available to the system at a given energy to the corresponding thermodynamic phase space. 
To this end, we note that all hadrons except pions are sufBciently heavy so that the Boltzmann approximation (55) 
for the particle spectra is applicable. In this case, the chemical potential only enters the overall normalization of the 
invariant momentum spectrum, but does not change its shape. For pions, on the other hand, the chemical potential 
is in general quite small [6-8] . It is therefore justified to use the following ansatz for the particle spectra 

1 dUn TO_L„COShy„ 

(56) 



m±n dm±ndyn " em±nCoshyn/T ^ I ' 

where c„ contains the information about the volume of the system, the degeneracy, and the chemical potential and is 
determined from the multiplicity of the nth particle, 

{rin) = Cn y" d^Pn ^jg^/j _^ ^ ■ (57) 

The rapidity spectrum drin/dyn is obtained from Eq. (56) by integrating over mj^„. 

In order to compute the invariant momentum spectrum (56), we still need to fix the temperature. We choose it 
such that the total thermal energy E is equal to the center-of-momentum energy ^/s available for particle production, 

V-S^E = ^[ d^PrE, = E^'- / '^'Pr , (58) 

^ J m±r dm±rdyr ^ J e^{Er/T) ±1 

where r labels the particle species. 



IV. COMPARISON WITH EXPERIMENTAL DATA 



In this section, we compute transverse mass and rapidity spectra from both the phase space as well as the thermal 
model and compare them to the experimental data. 



A. Completing particle multiplicities in central Pb+Pb collisions at 158 AGeV 

The calculation of the phase-space integral requires knowledge of the multiplicities of the various particle species, 
their masses, and the total energy of the reaction. Since the NA49 experiment does not measure all particle species 
produced in a nuclear collision, the multiplicity data [20-25] must be completed by using symmetry considerations, 
electric charge, baryon number and strangeness conservation, as well as empirical rules. 

We assume that the n particles in the final state are solely produced by strong interactions. Therefore, the 
multiplicities measured in the experiment have to be corrected for weak decays. This has already been done in the 
data of the NA49 collaboration [20]. We show the measured multiplicities for pions, kaons, A's, S's, and Q's in Table 
I. 
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4.45 


K+ 
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K- 
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n 
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0.62 


K% 


0.4977 


81.0 


Q. 


1.6725 


0.20 


A 


1.1157 


53.0 









TABLE I: Multiplicities of different particle species measured by the NA49 collaboration [20] in Pb+Pb collisions at 158 AGeV 
and corrected for weak decays. 



Based on Table I, we now reconstruct the multiplicities for particle species not directly measured by NA49. 

• 7r°: Due to (approximate) isospin symmetry in the reaction, we take 

{A = ^^^^^^ . (59) 

• K°,K'^: We assume that {K%) = {Kl) and {K°) = {K"). Since {K°) + = {K°) + {Kl), we have 

(if 0) = (i?0) = {Kl) . (60) 

• A, S+, E~ : An empirical rule found by Wroblewski [40, 41] states that 

(S+) + (E-)=0.6[(A) + (I]0)] . (61) 

Here, (A) + (E") corresponds to the number of measured A's as given in Table I. Note that the measured A 
multiplicity is not the true multiplicity, (A), because the E" decays to 100% into a A. One therefore has to 
reduce the measured A multiplicity by that of the E's. Equation (61) yields the combined multiplicity of E+ 
and E~ directly from the measured data. Due to isospin symmetry, we may assume that 

Now the multiplicities of E+ + E~ and E° are known. Combining Eqs. (61) and (62), we can solve for (A), 

(A) = I (EO) = I [(E+) + (E-)] = 0.7 [(A) + (EO)] . (63) 

The multiplicities for the corresponding antiparticles are obtained in the same way. 

• Due to isospin symmetry we have 

(SO) = (S-) , (SO) = (S+) . (64) 

• Nucleons: The net number of protons and neutrons follows from baryon number conservation, 

(j>) - (P) + (n) - (n) 

= N^- [(A) - (A) + (E+ + E-) - (E+ + S-) + (S°) - (5°) + (S") - (S+) + (17") - (17+)] 

= 276.96, (65) 

where we have used the number of wounded nucleons and the previously determined multiplicities. Assuming 
isospin symmetry, 

ip) = in) , iP) = in) , (66) 

we get (p) — (p) = (n) — {n) = 138.48. In order to obtain the individual multiplicities, we assume that the 

production of strange antiquarks relative to strange quarks is suppressed as compared to the production of 
nonstrange antiquarks relative to nonstrange quarks by a factor D = [42, 43]. This results in 

(p) (A) (S-) {n-) ^'"> 
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TABLE II: Completed list of multiplicities. 



D can be obtained from the measured multiplicities. We take the average value of / and |q-| / ' 
respectively. This yields D = 1.93. The antiproton-to-proton ratio is then ||y = 0.0454. 

The completed multiplicities are summarized in Table II. The average number of the secondaries thus obtained is 
2611. Note that we do not attempt to reconstruct resonances decaying strongly into pions and nucleons, such as 
the p and the A. In order to determine the corresponding multiplicities one would require additional assumptions, 
e.g. thermodynamical equilibrium between various particle species. As we aim to present an approach that is com- 
plementary to the thermal model, we refrain from doing so. Moreover, one would then have to treat the decay of 
these reconstructed resonances separately in order to compare with the measured spectra which already contain decay 
products. We leave this for a future improvement of our approach and, in this work, do not make a distinction between 
direct particle production and subsequent resonance decay. 



B. Transverse mass and rapidity spectra from phase space, o = 



In this section, we compute pion and kaon spectra from the phase-space model with a = and compare them to 
experimental NA49 data from the 5% most central Pb+Pb collisions at 158 AGeV. 

Figure 4 shows transverse mass spectra for 7r~, , and K~ mesons, calculated from the phase-space model as well 
as from a thermal model, together with the measured data. We observe that the phase-space model and the thermal 
model give quite similar results. The agreement becomes better for larger particle mass. This is not surprising, since 
the nonrelativistic limit of relativistic invariant phase space, J d^Pi/Ei — > (1/m,) / d^Pi, is identical to the phase 
space in the thermal model. 

However, both models disagree with the data. The slope is too small, i.e., the inverse slope parameter is too large. 

In the phase-space model, T* = 0.4935 GeV and T'^ = 0.4934 GeV. In the thermal model, the temperature parameter 
extracted from the energy constraint (58) is T = 0.374 GeV for both pions and kaons. Note that the inverse slope 
parameter T* is larger than the temperature T, T* /T ~ 1.32. In the ultrarelativistic limit, T*/T = |, cf. discussion 
above. In the nonrelativistic limit, i.e., for rrii ^ T, T* /T — s- 1, since then invariant phase space becomes identical 
to thermodynamic phase space. Pions and kaons are neither ultrarelativistic nor nonrelativistic, which explains the 
value T*/T ~ 1.32, i.e., between 1 and 3/2. In the range of rn± values shown in Fig. 4, both spectra appear to 
have a similar slope, but at larger m±, the thermal model spectrum decays faster than the one calculated from the 
phase-space model, as it should, since its inverse slope parameter T is smaller than T* by ^ 30%. 

Let us now consider the rapidity spectra. Figure 5 shows the measured spectrum of tt^ mesons in comparison to the 
phase-space model and the thermal model. Phase-space and thermal model give almost identical results. However, 
the calculated rapidity spectra are much narrower than the measured one. Similar results hold for mesons, see 
Fig. 6. 

How sensitive are the calculated spectra to variations in the total multiplicity? To answer this question, we 
performed two tests. First, we added 200 massless particles and, alternatively, 200 particles with mass 1 GeV, to 
the set of 2611 particles. For a given, sufficiently large y/s, this increases the derivative of ln/„_i with respect to 
^/s. Then, according to Eq. (48), the slope parameter 1/T* of the transverse mass spectra increases. This increase 
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FIG. 4: Transverse mass spectra for -k^ ,K^, and K~ calculated from the phase-space model with a = 0, as well as from a 
thermal model. The pion data is taken in the range < y < 0.2, and the kaon data for \y\ < 0.1. Experimental data is shown 
by the symbols. Here, as in the following figures, error bars include statistical errors (as given by NA49) and a systematical 
error of 5%. 



is independent of rapidity and the mass of the considered particle species. It amounts to 7.5% when adding massless 
particles and to 14% when adding particles with mass 1 GeV. The difference between adding massless and massive 
particles can be explained from Fig. 1: the increase of ln/„_i with ^/s is smaller for massless particles, thus, according 
to Eq. (48), the increase of the slope parameter 1/T* when adding particles is also smaller. An uncertainty of 200 
in the total multiplicity can therefore not explain the deviation between the measured and the calculated transverse 
mass spectra for pions and kaons. 

Second, we changed the total multiplicity by A; x 260 particles, where k = —1, 1, 2, preserving the ratio between the 
multiplicities of different particle species. Figure 7 shows the transverse mass spectra. The inverse slope decreases 
with increasing multiplicity, but even 560 additional particles are not sufficient to reproduce the measured spectra. In 
Fig. 8 we show the corresponding rapidity spectra for k = —1, 0, 1 together with the experimental data in a logarithmic 
plot. The distributions become narrower with increasing multiplicity. Thus, while the slope values of the transverse 
mass spectra approach the data for increasing multiplicity, the rapidity distributions deviate more strongly. It is 
therefore not possible to simultaneously reproduce transverse mass spectra and rapidity spectra simply by varying 
the total multiplicity. 



C. Transverse mass and rapidity spectra from phase space, a ^ 



We have seen that the transverse mass spectra calculated from the phase-space model with a = do not agree with 
the experimental data: their inverse slope parameter is too large. Similarly, the calculated rapidity spectra are too 
narrow. This indicates a substantial remnant of the initial longitudinal motion of target and projectile, which still 
dominates the phase space. Or in other words, complete isotropization has not occurred. 

The introduction of the factors fi{p±i), cf. Eq. (6), with a > 0, allows to suppress the transverse in favor of the 
longitudinal phase space. The parameter a regulates the strength of the suppression of transverse momenta and has 
to be extracted from the data. There are in principle two possibilities: (1) a is adjusted to reproduce the transverse 
mass spectra, or (2) a is adjusted to reproduce the rapidity spectra. 

In Figs. 9 and 10 we show the transverse mass and rapidity spectrum for tt" mesons, calculated from the phase- 
space model in comparison with experimental data. Choosing the value a = (0.2 GeV)~^, we can quite well reproduce 
the transverse mass spectrum. However, the rapidity spectrum is then wider than the data, which indicates that the 
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FIG. 5: Rapidity spectra for tt mesons, experimental data are shown by symbols, results from the phase-space model by the 
dashed line and from the thermal model by the solid line. 
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TABLE III: Optimal values for the damping parameter a. 



suppression of transverse phase space is slightly too large. Note that, at small transverse masses, the pion data lie 

above the calculated curve. This is quite naturally explained by the fact that our simple phase-space model does not 
account for the distortion of phase space due to resonance decays. Such resonance decays are expected to populate 
mainly the low-transverse mass region where we observe the diifcrence between the data and the calculation [15]. 

Figures 11 and 12 show the transverse mass and rapidity spectrum, now for a value a = (0.325 GeV)~^ which is 
adjusted to reproduce the rapidity spectrum. In this case, the calculated transverse mass spectrum differs from the 
data in that the inverse slope comes out too large. Transverse phase space is not sufficiently suppressed in this case. 
Comparing the set of Figs. 9, 10 with the set of Figs. 11, 12 suggests that there is a value for a between (0.2 GeV)~^ 
and (0.325 GeV)~^, for which a simultaneous fit of both transverse mass and rapidity spectrum, albeit of worse quality 
than the individual fits shown in Fig. 9 and Fig. 12, could be possible. Before returning to this question, we summarize 
in Table III the optimal fit values for the parameter a for 7r~, K^, and K~ mesons. In the first column, we quote the 
experimentally measured value of the slope parameter of the transverse mass spectrum, T*^^. In the second column 
we show the a-values optimized to fit the transverse mass spectra, and in the last column those optimized to fit the 
rapidity spectra. Note that, for the kaon spectra, the two a values thus obtained are quite close to each other. Indeed, 
either a fit to the transverse mass spectrum or to the rapidity spectrum will also yield a reasonable description of the 
respective other spectrum. 



Wc now ask the question whether there is a single value for a which gives a reasonable description of both the 
transverse mass and the rapidity spectrum, and not only for pions but also for kaons. This is indeed the case, as shown 
in Figs. 13, 14, and 15, for a = (0.256 GeV)~^. It is quite remarkable that, for a given set of particle multiplicities, a 
single parameter that accounts for longitudinal phase-space dominance is sufficient to yield a simultaneous fit of the 
transverse mass and rapidity spectra for both pions and kaons. The observed deviations can be quite easily explained. 



16 



80 




FIG. 6: The same as in Fig. 5 for K'^ mesons. 



In the case of the pions, resonance (e.g. p and A) decays populate the region of small transverse momenta and thus 
lead to an enhancement of the measured as compared to the calculated transverse mass spectrum. For the kaons, 
one observes that the measured transverse mass spectra exhibit a convex curvature indicative for collective flow of 
matter. This is a nontrivial dynamical effect not contained in our simple phase-space approach. Collective flow is 
more prominent for heavier particle species not considered here and will lead to similar deviations from the spectra 
calculated within our simple phase-space approach. Nevertheless, we expect that the quality of the fit can be further 
improved by introducing separate parameters for each particle species r. This, however, as well as a more refined 
treatment of resonance decays, is beyond the scope of the present work where we want to keep the treatment as simple 
as possible. 



V. SUMMARY AND OUTLOOK 



We have studied particle production in nuclear collisions within a phase-space model. Assuming a simple form 

for the matrix element that describes the production of n secondary particles, we have computed the phase-space 
integral within a statistical approach based on the Laplace transformation and the saddle-point method. Single- 
inclusive invariant momentum spectra can be then computed from the ratio of the phase-space integrals for n — 1 and 
n particles, respectively. 

We have applied this phase-space approach to multiparticle production to data for Pb-hPb collisions at 158 AGeV 
taken by the NA49 collaboration. In order to compute the phase-space integral, we require (a) the total available 
energy which is determined from the beam energy and the number of wounded nucleons (in our case, ^/s = 3131.3 
GeV) and (b) the multiplicities and the masses of the various particle species produced in the reaction. However, the 
NA49 experiment docs not mcasiire all species of strongly interacting particles. We reconstructed the multiplicities of 
the most important missing particle species using symmetry considerations, conservation laws, a simple quark model 
and an empirical rule. In this way we obtained 2611 secondary particles. By changing this value by - 10% and -|-20%, 
respectively, we tested the sensitivity of our results to the reconstruction procedure and found that this does not affect 
our conclusions. 

Prom the phase-space integral we then computed transverse mass and rapidity spectra for tt", K^, and K~ mesons 
and compared them with the measured spectra. Our simple model predicts transverse mass spectra of exponential 
shape, with an inverse slope parameter that is (almost) independent of the particle mass. 

We also computed transverse mass and rapidity spectra assuming the system to be in thermal equilibrium at a 
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temperature which is determined from the total available energy y/s = 3131.3 GcV. In the thermal model, transverse 
mass spectra are not pure exponentials, but are of slightly convex shape. At midrapidity, the inverse slope parameter 
of the transverse mass spectra is identical to the temperature, T, and smaller than the inverse slope parameter, T*, 
arising in the phase-space model. This difference can be traced to the different definitions of phase space in the 
two models, invariant phase space d^pi/Ei for the phase-space model, and non-invariant phase space (Ppi for the 
thermal model. For ultrarelativistic particles, the difference in the slope parameters can be computed analytically, 
T*/T = 3/2. Due to the convex shape of the transverse mass spectra in the thermal model, in the range of transverse 
masses < m± — m < 1 GeV, the transverse mass spectra in both models are surprisingly similar. The difference in 
the slope parameters becomes visible only at larger transverse masses. The rapidity spectra of the thermal and the 
phase-space model are also very similar. 

If we assume a constant, i.e., isotropic matrix element for multiparticle production, the transverse mass and rapidity 
spectra calculated from the phase-space model deviate substantially from the experimentally measured spectra. The 
same holds for the thermal model. This deviation is due to the fact that there is still a remnant of the initial, 
longitudinal motion of the nuclei. 

There are two possibilities to obtain agreement between models and data. The first possibility is to retain the 
assumption of an isotropic phase-space distribution locally in coordinate space. One then employs a dynamical model 
for the evolution of matter, such as fluid dynamics, which gives rise to collective transverse as well as longitudinal 
flow of matter. For a given equation of state, one can then tune the initial as well as the freeze-out conditions in order 
to reproduce the data and thus account for longitudinal phase-space dominance. 

In this paper, we pursue the second possibility. We relax the assumption of isotropic phase space and take longitu- 
dinal phase-space dominance into account by a simple ansatz for the matrix element which exponentially suppresses 
transverse momenta. This has the effect that the transverse mass spectra decrease faster and the rapidity spectra 
become wider. The amount of suppression is determined by a parameter a which is fitted to the experimental spectra. 
The best fits were obtained by choosing different parameters a for each type of particle (77+ , , K^) and fitting 
transverse mass spectra and rapidity spectra separately. Surprisingly, a rather good fit of transverse mass and ra- 
pidity spectra of all particles considered can also be obtained by a single value, a = (0.256 GeV)"'^. Note that this 
second approach does not require the assumption of thermal equilibrium. Observed deviations between calculated 
and measured spectra can be attributed to nontrivial dynamical phenomena, such as resonance decays and collective 
flow. 

A possible extension of our present work is to compare our phase-space model with data for heavier particles, such 
as nucleons and hyperons. It is known that the transverse mass spectra of particles heavier than pions and kaons have 
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FIG. 8: Rapidity spectrum for tt mesons for different total multiplicities n. 



a substantially larger inverse slope parameter (but essentially of similar magnitude for all of them). This is usually 
attributed to the collective radial flow of matter [44]. In order to incorporate this into our approach, we would have 
to introduce one additional parameter a for the heavier particle species. The other way to account for different slope 
parameters for the various particle species would be to use a more physical ansatz for, or even a detailed calculation 
of, the matrix element T(pi) [17]. Another possible extension is to incorporate the distortion of the phase-space of 
light particles, such as pions, by the decay of heavier resonances [13]. 
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FIG. 9: Transverse mass spectrum for n mesons in the range < y < 0.2. The sohd hne shows the spectrum calculated from 
the phase-space model for a = (0.2GeV)~^. 
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FIG. 10: The same as in Fig. 9, but for the vr rapidity spectrum. 
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FIG. 14: Rapidity spectrum for tt mesons. The solid line shows the spectrum calculated from the phase-space model with 
a = (0.256 GeV)-\ 



